Standing and Moving Gap Solitons in Resonantly Absorbing Gratings 
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We present hitherto unknown forms of soliton dynamics in the forbidden frequency gap of a Bragg 
reflector, modified by periodic layers of near-resonant two- level systems (TLS). Remarkably, even 
extremely low TLS densities create an allowed band within the forbidden gap. This spectrum gives 
rise, for any Bragg reflectivity, to a vast family of stable gap solitons, both standing and moving, 
having a unique analytic form, an arbitrary pulse area, and inelastic collision properties. These 
findings suggest new possibilities of transmission control, noise filtering or "dynamical cavities" 
(self-traps) for both weak and strong signal pulses. 
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The study of light-matter interactions in periodic di- 
electric structures has developed into a vast research 
area. At the heart of this area is the interplay between 
the resonant reflections induced by the Bragg reflector, 
giving rise to photonic band gaps, and their dynamical 
modifications due to nonlinear light-matter interactions. 
The pulsed mode of propagation in such structures ex- 
hibits a variety of unique fundamentally and technolog- 
ically interesting regimes: nonlinear filtering, switching, 
and distributed-feedback amplification Of particu- 
lar interest are gap solitons (GS), i.e., moving or stand- 
ing (quiescent) self-localized pulses, whose spectrum is 
centered in a gap induced by the grating. GS in Kerr- 
nonlincar Bragg reflectors have been extensively analyzed 
Q and experimentally observed ||. Recently, GS have 
also been predicted in Bragg-reflecting second-harmonic 
generating media B. 

This work is dedicated to a different mechanism sup- 
porting GS in periodic media, which is based on near- 
resonant field-atom interactions. The first step in this 
direction has been made in Ref . 0] , where an exact mov- 
ing GS solution has been found in a periodic structure 
composed of thin layers of resonant two-level systems 
(TLS) separated by half-wavelength non-absorbing di- 
electric layers, i.e., a resonantly-absorbing Bragg reflector 
(RABR) . In the soliton solution obtained in Ref. B , the 
combined area of the forward- and backward-propagating 
pulses is 2n, characteristic of self-induced transparency 
(SIT) solitons in uniform media ||. The existence of 
this soliton stems from the cooperative resonant atomic 
polarizability, which compensates for the periodic mod- 
ulation of the linear polarizability in the Bragg reflector 
1^,0. The analysis presented in Ref. || leaves several 
open questions of fundamental and applied importance: 
(z) Can one overcome the basic restriction implicit in this 
solution, namely, that the cooperative length over which 
a SIT pulse is formed, must be shorter than the Bragg 
reflection length? If this restriction is essential, then the 



soliton would only exist in weakly-reflecting Bragg struc- 
tures, which can hardly serve as efficient filters that block 
pulses other than GS. (ii) Are GS admitted in a RABR 
for weak pulses whose area is less than 2tt? (iii) Is there a 
quiescent counterpart to the moving GS, which would im- 
ply complete dynamic confinement of light in the RABR? 
(iv) What is the result of collisions between moving GS 
in this system? 

In this work, we present answers to the above ques- 
tions: (a) The RABR supports a vast family of GS char- 
acterized by two parameters, the soliton amplitude and 
velocity (analogously to GS in Kerr-nonlinear gratings 
0), which exists for any ratio of reflection to cooperative 
lengths, (b) It includes a subfamily of quiescent or zero- 
velocity (ZV) solitons with an arbitrary pulse area, whose 
exact analytical form represents an essentially novel type 
of soliton solutions in nonlinear optics, (c) Moving GS 
are found analytically in the small- amplitude limit, where 
they reduce to the nonlinear-Schrodinger (NLS) solitons, 
and also numerically, as deviations from the exact sub- 
family found in ^ or from the exact ZV solitons. Our 
simulations indicate that both ZV and moving GS are 
stable, (d) We simulate collisions between moving GS, 
demonstrating that they may be either weakly or strongly 
inelastic. These findings reveal a basically new phe- 
nomenon: the multiple reflections in a grating can effec- 
tively change the pulse coupling to the TLS, so that soli- 
ton transmission is not restricted to pulses of 2tt area, as 
in ordinary SIT, and occurs (at an appropriate frequency) 
for any ratio of reflection to cooperative lengths. 

Our starting point are the equations for the sum E + 
of the amplitudes of the forward- and backward- propa- 
gating electromagnetic waves, polarization P, and pop- 
ulation inversion w, derived in Ref. J^] from the coupled 
Maxwell-Bloch equations for the bidirectional propaga- 
tion in RABR: 
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where r and £ are the normalized time and propaga- 
tion distance, 8 is the effective detuning of the field from 
the atomic resonance. The key parameter 77 = l c /l r is 
the ratio of the Arrechi-Courtens cooperativity length 
(L = 2c/uj p , where uj p is an effective plasma frequency 
p|) to the reflection length (l r = Ade^j V Ae, where d is 
the period and Ae the variation of the dielectric index of 
the periodic structure with average dielectric index eo 

f. The crucial assumption in the derivation of Eqs 
is that the resonant absorbers are confined to thin lay- 
ers, periodically inserted between passive thick dielectric 
layers. The system (0)-(|) can be simplified: substitut- 
ing into Eq. (§) 

the expression for E+ following from (H), 
one obtains an equation that can be explicitly integrated 
to yield a simple algebraic relation 
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w = ±y 1 
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that should be further inserted into Eq. ([!]). 

First, we linearize Eqs. in order to obtain a 

crucially important characteristic of the model, its dis- 
persion relation, by fixing w = — 1 and substituting into 
the linearized equations ~ exp (ik( — iur) . The 

resulting dispersion relation: 

(w - S) [u, 2 -k 2 -(2 + v 2 )] + 2 (77 - 6) = 0, (5) 

is displayed in Fig. [I]. The frequencies corresponding to 
k = arc w = r] and w = — \ (77 — S) ± ^2+ j(r] + 5) 2 , 

while at k 2 — > 00 the asymptotic expressions for differ- 
ent branches of the dispersion relation are to = ±fc and 
uj = S + 2 (r/ — S) k~ 2 . Thus, the linearized spectrum al- 
ways splits into two gaps, separated by an allowed band, 

except for the special case, 77 
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when the upper gap closes down. The upper and lower 
band edges are those of the periodic structure, shifted by 
the induced TLS polarization in the limit of a strong re- 
flection. They approach the SIT spectral gap for forward- 
and backward-propagating waves R] in the limit of weak 
reflection. The allowed middle band corresponds to a 
collective atomic polarization excitation. 

Inside these gaps, ZV solitons are sought in the form 
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with real functions cr(C) and (?(£), X being the frequency 
detuning from the gap center. Using Eqs. ^ and (Q), 
we can express the variables w and q in terms of a: 
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The remaining equation for <t(£) is 
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where the ± signs correspond to those in Eqs. (^). 

Equation (||) is the Newton equation of motion for a 
particle with the coordinate a in the potential F(a). It 
gives rise to a soliton-like solution (|], provided that the 
upper sign is chosen in Eqs. (0), and the frequency x 
satisfies the conditions |%| > 77, \(x + v){x ~ $)\ < 2. 
These inequalities can be solved to yield, in an explicit 
form, two frequency intervals in which one has the ZV 
gap solitons: at x > 0, it is 



V<X<\J\{v + 5) 2 + 2- l -{r 1 -5), (9) 

which only exists for small and moderate 77 (the weak 
Bragg-reflectivity limit) 



r)<~6 + Jl + j<5 2 , 
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and another interval at x < 0, which exists for any re- 
flectivity 



V < -x < 
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On comparing these expressions with the spectrum 
shown in Fig. we conclude that part of the lower gap is 
always empty from solitons, while the upper gap is com- 
pletely filled with ZV solitons in the weak-reflectivity case 
( |l0| ) (Fig. 0(a)), and completely empty in the strong- 
reflectivity case (JTTJ) (Fig. 0(b)). It is relevant to mention 
that a partly empty gap has recently been found in the 
model combining Bragg reflection and second harmonic 
generation 

Inside the frequency intervals (||) and (|ll|) , Eq. (JsJ) can 
be integrated by means of the substitution 
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This yields the soliton shape in an implicit form, i.e., £ 
vs. p: 



a(0, P(Cr)=ie-^q(0, (6) | C | = 2 
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where p\ = 1 — \ \ (x + v) (x ~ $) I (note that p\ is positive 
under the above conditions (^)-(0)). It can be checked 
that this ZV gap soliton is always smgZe-humped. Its 
amplitude can be found from Eq. (p"3|), 
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The most drastic difference of these new solitons from 
the well-known SIT pulses || is that the area of the ZV 
soliton is not restricted to 2tt, but, instead, may take 
an arbitrary value. As mentioned above, this basic new 
result shows that the Bragg reflector can enhance (by 
multiple reflections) the field coupling to the TLS, so as 
to make the pulse area effectively 2n. In the limit of 
the small-amplitude and small-area solitons, p$ -C 1, Eq. 
( |l3| ) can be easily inverted, the ZV soliton becoming a 
broad sech-likc pulse: 



% - %o sech 
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(nl) remains finite, but the soliton width diverges 
(pig. |(b)). Thus, although the ZV soliton 



In the opposite limit, 1 — Pq — > 0, i.e., for vanishingly 
small \x + tne soliton's amplitude becomes very 
large, and further analysis reveals that, in this case, the 
soliton is characterized by a broad central part with a 
width ~ (l — p^) 1//2 (Fig. ^](a)). Another special limit 
is X~ V ~ * 0- It can be checked that in this limit, the am- 
plitude 
as |x-i7| 

has a single hump, its shape is, in general, strongly dif- 
ferent from that of the traditional nonlinear-Schrodinger 
(NLS) sech pulse. 

The stability of the ZV gap solitons was tested numer- 
ically, by means of direct simulations of the full system 
(Q) - (H), the initial condition taken as the exact soliton 
with a small perturbation added to it. Running the simu- 
lations at randomly chosen values of the parameters, we 
have always found the ZV GS to be apparently stable. 
However, possibility of their dynamical and structural 
instability need to be further investigated, as has been 
done in the case of GS in a Kerr-nonlinear fiber with a 
grating |,@. 

Although the system of Eqs. (0) - (||) is not explic- 
itly Galilean- or Lorentz-invariant, translational invari- 
ance is expected on physical grounds. Hence, a full fam- 
ily of soliton solutions should have velocity as one of 
its parameters. This can be explicitly demonstrated in 
the limit of the small-amplitude large-width solitons [cf. 
Eq. (p"5|)]. We search for the corresponding solutions 



in the form £+(£, t) = <r(£, r) exp (— iw r) , P(£, r) 



i(j(C, t) exp (— iwor) (cf. Eqs. (g)), where u>o is the fre- 
quency corresponding to k = on any of the three 
branches of the dispersion relation (||) (see Fig. |l|), and 
the functions a((, r) and q((, r) are assumed to be slowly 
varying in comparison with exp (— itu^r). Under these as- 
sumptions we arrive at the following asymptotic equation 
for it(C,t): 

.lo q {lo q - S) 2 - i] + 5 d d 2 
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Since this equation is of the NLS form, it has the full 
two-parameter family of soliton solutions, including the 
moving ones S. 

In order to check the existence and stability of the 
moving solitons numerically, we simulated Eqs. ([!]) - (|^) 
for an initial configuration in the form of the ZV soliton 
multiplied by exp(ip^) with some wavenumberp, in order 
to "push" the soliton. The results demonstrate that, at 
sufficiently small p, the "push" indeed produces a moving 
stable soliton (Fig. 0(c)). However, if p is large enough, 
the multiplication by exp(ip£) turns out to be a more 
violent perturbation, splitting the initial pulse into two 
solitons, one quiescent and one moving (Fig. ||(d)). 

A one-parameter subfamily of moving GS was found 
in an exact form in Ref. [pi : 



Aq exp [i (a( — At)] sech 
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where A is the detuning from the gap center, and the 
squared am- 
plitude A 2 = (1 - u 2 )' 2 \&u 2 (1 - u 2 ) - r, 2 (1 + u 2 ) 2 } is 

expressed in terms of the velocity u (normalized to c). 
The values of u are restricted by the condition Aq > 0, 
which, in particular, forbids u — 0. The present analy- 
sis strongly suggests, but does not rigorously prove, that 
the subfamily (|l7]) belongs to a far more general two- 
parameter family, whose other particular representatives 
are the exact ZV solitons ( |l3| ) and the approximate small- 
amplitude solitons determined by Eq. (fi6|). 

An issue of obvious interest is collisions between GS 
moving at different velocities. In the asymptotic small- 
amplitude limit reducing to the NLS equation (|l6|), the 
collision must be elastic. To get a more general insight, 
we simulated collisions between two solitons given by 
(|l7j). The conclusion is that the collision is always in- 
elastic, directly attesting to the nonintegrability of the 
model. Typical results are displayed in Fig. [3| which 
demonstrates that the inelasticity may be strong, de- 
pending on the parameters. 

The present findings can be demonstrated in peri- 
odically etched dielectric structures containing cither 
strained quantum wells or gas as the active TLS. For ex- 
ample, we can use HF gas whose active dipole transition 
at the wavelength A = 84 /zm is resonant with a pho- 
tonic band edge. A 4.5 mTorr gas pressure corresponds 
to cooperative length l c = 114 cm and inhomogeneous 
dephasing length cT 2 * ~ 10 4 cm. In a structure with 
Bragg reflection length l r = 100A <C l c , we then have 
ri = -S= -135.7, which allow for GS detuned by 2 x 10 3 
cm -1 from both band edges with ~ 2 ns width and 4x 10 9 
s _1 Rabi frequency. The suggested structures can be 
used to realize any kind of GS: standing (ZV) GS require 
initial localized excitation of the near-resonant medium 
by counterpropagating laser beams, whereas moving GS 
can be launched by a single near-resonant laser beam 
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propagating along the structure. 

In conclusion, we have demonstrated that a periodic 
array of near-resonant two-level systems (TLS) combined 
with a Bragg grating gives rise, for any Bragg reflectiv- 
ity, to a vast family of stable gap solitons, both standing 
and moving, having a unique analytic form, an arbitrary 
pulse area (Fig. ||), and inelastic collision properties (Fig. 
||). Remarkably, even extremely low TLS densities cre- 
ate an allowed band within the forbidden gap (Fig. [I]). 
These findings reveal hitherto unknown forms of soliton 
dynamics in a nonintegrable, strongly nonlinear system. 
Their origin is the surprisingly rich interplay between 
multiple reflections and cooperative, near-resonant field- 
matter interaction, which removes the pulse-area restric- 
tions of ordinary self-induced transparency. These find- 
ings can lead to the realization of novel, highly advan- 
tageous filters, which can stably transmit selected signal 
frequencies through their spectral gaps, while effectively 
blocking others, with no restriction on the signal pulse 
area. Alternatively, they can be used to spatially confine 
(self-trap) light in certain frequency bands, thus creating 
"dynamical cavities" . 
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FIG. 1. The spectrum of the linearized RABR model 
([l|)-(|^). The parts of the upper and lower gaps filled with 
the solitons are shaded: (a) the weak-reflectivity case when 
the condition @ holds (jj = 0.6, 8 = 0); (b) the opposite 
case (rj = 2, 6 = 0). 
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FIG. 2. Zero- velocity solitons |£+(C)| 2 (a) S = 0, n = 0.9, 
X = —0.901 (divergent width and amplitude); (b) idem, but 
for x = 0.901 (divergent width and finite amplitude), (c) 
Pulse(s) obtained as a result of "pushing" of zero-velocity 
solitons (dashed lines), by the initial multiplier exp(— ipC,) af- 
ter a sufficiently long evolution (r = 400) (solid lines). 5 = 0, 
V — 4, X — —4.4, and p = 0.1. (d) idem, but for p = 0.5. 
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FIG. 3. Typical example of inelastic collisions between the 
solitons (jlT]) at 5 — and r\ = 0.5, with the velocities (nor- 
malized to c) u\ = 0.6, U2 = —0.75. 
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